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1. Introduction 

The First ISCCP Regional Experiment (FIRE) Cirrus Intensive Field 
Observation (IFO) program measured cirrus cloud properties with a variety of 
instruments from the surface, aircraft, and satellites. Surface and 
aircraft observations provide small-scale point and line measurements of 
different micro- and macrophysical properties of advecting and evolving 
cloud systems. Satellite radiance data may be used to measure the areal 
variations of the bulk cloud characteristics over meso- and large - scales . 
Ideally, the detailed cloud properties derived from the small-scale 
measurements should be tied to the bulk cloud properties typically derived 
from the satellite data. Full linkage of these data sets for a 
comprehensive description of a given cloud field, one of the goals of FIRE, 
should lead to significant progress in understanding, measuring, and 
modeling cirrus cloud systems. In this paper, the relationships derived 
from intercomparisons of lidar and satellite data by Minnis et al . (1989a) 
are exploited in a mesoscale analysis of the satellite data taken over 
Wisconsin during the Cirrus IFO Case Study. 


2. Data and Methodology 

Daytime half-hourly, 1-km visible (VIS, 0.65 pm) and 4x8-km infrared (IR, 
11.5 pm) GOES data were used to construct bispectral histograms on a 0.5° 
latitude - longitude grid between 42°N and 47 N and 87 W and 92 W during 27-28 
October 1986. The VIS data are stored as counts between 0 and 63 which were 
converted to radiance and reflectance, p. The IR data are stored as counts 
between 0 and 255 which are converted to equivalent blackbody temperature, 

T. Soundings taken every 6 hours over Green Bay, Wisconsin were used to 
determine the vertical variation of temperature. Linear interpolation 
between the soundings was used to estimate the temperatures at each half 
hour . 

The parameters of interest here are cloud amount, C; cloud- top 
temperature, T ; cloud- center temperature, T c ; cloud- top height, z t ; cloud- 

center height, z c ; cloud thickness, h; cloud emittance, e; VIS optical 

depth, r , and clear-sky temperature, T cs - The relationships between these 

parameters during the IFO were explored by Minnis et al . (1989a, b) using 
combined satellite-lidar data. A simple physical model is used here to 
relate the observed reflectance to cloud optical depth: 
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Cloud albedo is q - a (r p 0 ) ; the diffuse cloud albedo is a, - a At )• 

ccv ci a v v 

and the cloud anisotropic reflectance factor is x q • The values for the 

albedos are determined in the same manner as Minnis et al . (1989b) using the 
results of Takano and Liou (1989). Values for x c are taken from the 

results of Minnis et al . (1989a). Ozone transmittance above the cloud is 
T a’ T c and T u are the transmi - ttances of the cloud to direct downward and 

upward VIS radiation, respectively; p is the clear-sky reflectance- a 

s ’ sd 

is the diffuse clear-sky albedo; and p 0 is the cosine of the solar zenith 
angle. The infrared optical depth is t & - 2.17r . The effective infrared 
beam emittance is e - 1 - exp(-T g / p) , where p is the cosine of the 
viewing zenith angle. Assuming a pixel is completely cloud-filled, 

T c - 1 { [ B(T) - (1 - OB(T cs )] / O. (2) 

Details for computing T z z and h are given in Minnis et al . 
(1989b). 

The basic approach to deriving these parameters from a histogram is 
outlined below. 

1) Consider all pixels with T < T - 3K as cloudy. 

2) Assign pixel to altitude level: low, middle, or high (Fig. 1). 

3) Compute average reflectance for each temperature in a layer. 

4) Compute for each observed temperature within the layer. 

5) Compute T c , T , and e for each observed T in the layer. 

6) Compute average T fi , T t> t^, and e for each layer. 

7) Compute z^, z t> and h for each layer from results of (6). 

8) Compute mean values for each variable by combining layer reults. 

The functions and describe the VIS-IR thresholds between low and 

middle and between middle and high cloud, respectively. 

3. Results and. Discussion 

This approach was applied to the lidar-satellite data of Minnis et al . 
(1989a) to determine the errors in the resulting cloud parameters. A 
comparison of selected parameters for total cloudiness is shown in Fig. 2 
for Fort McCoy data. Figure 3 shows the comparison of the lidar cloud- top 
heights and the satellite-derived high cloud- top altitudes. The comparison 
results are summarized in Table 1. The satellite -derived total cloud-center 
heights agree well with the lidar, on average, while the corresponding 
cloud- top altitudes are underestimated by 0.5 km. If only high clouds (as 
detected by the satellite) are considered, the satellite results agree with 

the lidar values of z , but underestimate the values of z These 

C c ' 

differences may be attributed to several effects. Some midlevel cloudiness 
may have been within the satellite field of view and still not have advected 
over the lidar sites. Partially cloud- filled pixels may have decreased the 
altitudes of some pixels or the values of x may have been inadequate at 

some hours. 


346 


I 



VISIBLE COUNT, D 


Fig- 


Examples of the derived values of high and midlevel cloud fraction, z^ , 

and £ are shown in Figs. 4-7 for 15, 18, and 21 UTC during October 28, 
1986. At 15 UTC, most of the clouds are above 6 km (Fig. 4). Midlevel 
cloudiness is confined to the southern portion of the region (Fig. 5). 
Average values of z range between 7 and 9 km (Fig. 6). The thicker 

clouds, as determined by the emittances, are primarily in the northern half 
of the IFO box (Fig. 7). By midday (18 UTC), an area of clearing has moved 
into central Wisconsin. Cloud-center heights have increased to 11 km in the 
east while some scattered high and midlevel clouds follow the clearing line. 
At 21 UTC, the average cloud heights vary from 7 to 9 km again with 
significant midlevel cloudiness in the northern half of the box. Some 
mesoscale structure is evident in the emittance fields. 

4 . Concluding Remarks 

The cloud parameters derived with the empirical technique described 
represent the most accurate, detailed areal cloud properties available for 
the Case Study. Additional research using other IFO results with 
theoretical calculations is needed to generalize the technique used here. 
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1. Schematic diagram of histogram analysis 

Table 
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0.94 

* t (km) 
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1. Differences between cloud parameter 
values derived with Ildar-satellite 
data set and satellite data only. 
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Fig. 3. Comparison of satellite-derived high cloud-top heights 
with lidar data for October 27-28, 1986. 
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